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FINITE DIFFERENCE ANALYSIS OF ROTORDYNAMIC SEAL 
COEFFICIENTS FOR AN ECCENTRIC SHAFT POSITION 


R. Nordmann and F.J. Oietzen 
Department of Mechanical Engineering 
University of Kaiserslautern 
Kaiserslautern, Federal Republic of Germany 


The dynamic coefficients of seals are calculated for shaft movements around an 
eccentric position. The turbulent flow is described by the Navier-Stokes equations 
in connection with a turbulence model. The equations are solved by a 
finite-difference procedure. 


INTRODUCTION 

To model the dynamic behaviour of turbopumps properly it is very important to 
consider the fluid forces which are developed in the seals. This has been clearly 
demonstrated by some authors, like for example Diewald (1987). The fluid forces are 
normally described by the following equation. 
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But this equation is only valid for a shaft moving around the center of the seal, 
which very seldom occurs in reality. In most machines the shaft will orbit around an 
eccentric position, so that the fluid forces must be described by 


fF 1 


K 

k 

r z 


D 

d i r z 

1 M o 1 f z I 

z 


zz 

zy 



zz 

zy • 

zz 

F 

= 

-k 

K 

1 Y 

+ 

-d 

D J 1 Y 

+ o M J 1 Y I 

y J 


y Z 

yy J 

L J 


yz 

yy J L - 

L yy J L J 


( 2 ) 


The dynamic coefficients in such a case have been investigated by Jenssen (1970), 
Allaire et al. (1976) and recently in an excellent paper by 
Nelson and Nguyen (1987). To model the turbulent flow, all these have used so 
called "Bulk-Flow Theories" in which the shear stress at the wall is described as a 
function of the average fluid velocity relative to the wall. Because some authors 
affirm that in the case of great eccentricities recirculation in circumferential 
direction occurs (which can’t be described by a bulk-flow model) and strongly 
effects the dynamic coefficients, we extended the theory of Dietzen and 
Nordmann (1987) and Nordmann (1987) to investigate the flow and the coefficients in 
the case of an eccentric shaft. In that theory the Navier Stokes equations in 
connection with the k-e turbulence model were used to calculate the dynamic 
coefficients of incompressible and compressible seals for a shaft motion around the 
centric position. 
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GOVERNING EQUATIONS 


In a turbulent flow, turbulent stresses occur which are often modeled like laminar 
stresses by introducing a turbulent viscosity. The turbulent and the laminar 
viscosity are then added to an effective viscosity. 

fx e = f^i + U t (3) 

The turbulent viscosity must be described by a turbulence model. We use the k-e 
model, but also much simpler mixing-length models will be appropriate for a 
straight seal. The turbulent viscosity is given by 
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So we have the Navier-Stokes equations, the continuity equation and the equations 
of the k-e model to describe the turbulent flow in a seal. These equations have the 
following form. 
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(The constants of the k-e model are given in appendix A) 
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PERTURBATION ANALYSIS 
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Fig. 1 Geometry of the eccentric shaft 


To describe the flow if the shaft is moving around an eccentric position we follow 
the procedure of Dietzen and Nordmann (1987) and use a similar transformation, 
(see Fig. 1) 


r (0) - r 

1 - 'aW ' -tfeO) — c ° (e) (6) 

But now r and Co are functions of 0. This is not so if the shaft moves around the 
center. 6 is the seal clearance, varying with angle 0 and time t. 

The radius r can be described by the following equation: 
a 


r ft (0) = - ex 2 sin 2 (0) - ex cos(0) (7) 

As result of this transformation the shaft orbiting around an eccentric position is 
transformed to the stationary eccentric position. 

If we introduce this transformation into our equations, we must obey the following 
relations. 
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To calculate the rotordynamic coefficients we assume that the shaft moves around 
the eccentric position on small orbits, so that we are allowed to introduce a 
perturbation analysis. 

6 = Co - eh t u = u Q + eu t v = v Q + ev t (9) 

w = w 0 + eWi p = Po + ep 4 

If we introduce these expressions and the coordinate transformation in our 
governing equation, neglecting terms with power of e greater than 1 and separating 
the parts with and without e we will get a set of zeroth order and first order 
equations. The zeroth order equations describe the stationary flow for the 
eccentric shaft, the first order equations the perturbation of this flow, if the 
shaft moves around the stationary position. 

We assume that the shaft moves with the precession frequency Q on a circular orbit 
with radius r 0 around the eccentric position. So the change of the clearance is 
given by 


eht = r 0 (cosOt cos0 + sinflt sin0) (10) 

Because this change is periodic in time we introduce also periodic functions for 
the flow variables. 

Uj = u Jc cosQt + u lg sinfit v t = v 1{; cosOt + v lg sinOt 

w t = w t cosOt + Wi sinOt Pt = Pj cosOt + p t sinOt 

c s OS 

By separating now in the first order equations the terms with cosOt and sinOt we 
obtain two real equations for every first order equation. These equations are 
then arranged in a new form by introducing complex variables. 
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Finally supplementary to our real zeroth order equations we have a set of complex 
first order equations. These equations have the following form. 
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Table 2: Source terms of zeroth and first order equations. 


Only the first order continuity equation to calculate p d has a slightly different 
form. 

* p5e <p " l) - D Po < 14 > 

You get S , S , S , S. , S , if you replace in table 1 in the corresponding 
Uq Vq Wq Kq 

terms r by rj and u,v,w,p,k,e by u 0 ,v 0 .Wo.Po>l<o.& 0 ' The terms C . C , C result 

u i v i w i 

from the perturbation of the convective terms in the Navier Stokes equations. 

D , D . D , D are constants resulting from the coordinate transformation, 
u 0 Vq Wo Po 

/n A A A 

which are not functions of u t . v lt w t , p* . The terms with H represent the time 
dependent parts. 
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Because we assume that the viscosity remains constant for the small motions, we 

A Aw 

do not need a k t and e t equation. 

(D . D , D . D . C , C . C are given in appendix B) 

Uq Vo W 0 Po Vj w t 


BOUNDARY CONDITIONS 

The zeroth order boundary conditions are 
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u^ is the average axial entrance velocity for every plane with 0 constant. Pjj es 
the reservoir pressure and u the rotational frequency of the shaft. 
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The first order boundary conditions are 


stator : 

u lg = (0.0) v lg = 

= (o.o) ; ls = (o.o) 
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THE FINITE DIFFERENCE METHOD 

For solving these equations a finite-difference procedure is used which is based on 
the method published by Gosman and Pun (1974). The seal is discretized by a grid 
(Fig. 2) and the variables are calculated at the nodes. The velocities u 0 ,v 0 . w 0 

AAA 

(Ui.Vt.Wt) are determined at points which lie between the nodes where the variables 

A 

Po.ko.^o (Pi) are calculated (Fig. 3) 
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Fig. 2 3-dimensional grid Fig. 3 Velocities and pressure 

in the seal. in a staggered grid. 

To calculate the flow we proceed as follows: 

1. We start our procedure with guessed values for all variables. 

AAA 

2. First the velocities u 0 ,v 0 ,w 0 (u^v^Wj) are calculated. 

3. Then the velocities and the pressure are corrected to satisfy 
the continuity equation. 

To do that we use a modified version of the ’PISO* procedure 
of Benodekar et al.(1985). 

4. After this k and e are calculated (only for the zeroth order solution). 

We repeat step 2 to 4 until we reach a convergent solution. First we solve the 
zeroth order equations and then the first order equations. 

Of course we need a 3-dimensional finite-difference method to calculate the flow in 
the case of an eccentric shaft, while a two-dimensional method is sufficient for 
movements around the centric position. 


DETERMINATION OF THE DYNAMIC COEFFICIENTS 

A 

By integration the pressure p t around the shaft we get forces in z and y directions. 

* • • ♦ • • 

Then we introduce z, y, z, y, z, y from our circular orbit into equation (2). 
This gives us the following equations: 


275 



L 2 r 


JI 
0 0 
L 2t 

JJ 
0 0 
L 2tt 


Pi- 


Pi. 


JJ 

0 0 

L 2n 

JJ 


Pi. 


Pi. 


cos0 R^ d0 dx 

cos0 R^ d0 dx 

sin0 R. m d0 dx 
lN 

sin0 R^ d0 dx 


- C °N < K zz * V - M zz° ! ) 
’ C °N < k zy - D zz°> 

= C OM (-k + D fi) 

°N v yz yy * 

= c OM (k + d n - m a 2 ) 

°N v yy yz yy J 


0 0 

If we calculate the forces for several precession frequencies Q of the shaft, we 
can obtain the coefficients by a ’Least-Square-Fit * . 


RESULTS 

We compare our theory with the model of Nelson and Nguyen (1987) and some 
experimental results of Falco et al. (1986) which also have been published in the 
paper of Nelson and Nguyen (1987). 

In Fig. 4-10 dynamic coefficients are calculated as a function of the eccentricity. 
We compare our results with Nelson and Nguyen’s theory and the stiffness 
coefficients also with experimental and theoretical results of Falco et al. 

The seal data are 

length : L = 40.0 mm pressure drop : 1.0 Mpa 

shaft radius : R^ = 80.0 mm shaft speed : 4000 RPM 

nominal clearance : C Q ^ = 0.36 mm 

density : p = 1000 kg/m 3 

preswirl ratio : w (O,0 )/Rjj^ = 0.3 

viscosity : = 1.0*10“ 3 Ns/m 2 

entrance lost-coefficient : f = 0.5 


CONCLUSIONS 

We have shown that it is possible to calculate the dynamic coefficients of 
eccentric seals by a finite difference method based on the Navier Stokes equations. 
This method can also be extended to calculate the coefficients of eccentric gas 
seals by following the procedure of Nordmann (1987) and to calculate the dynamic 
coefficients of bearings by neglecting the turbulence model. 
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Fig. 5 Direct stiffness versus eccentricity. GX/ C Q 
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Fig. 6 Cross coupled stiffness k versus eccentricity 
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Fig. 7 Cross coupled stiffness k versus eccentricity 
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NOMENCLATURE 
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Forces on the shaft in z and y direction 

direct and cross-coupling stiffness 
direct and cross-coupled damping 
direct and cross-coupling inertia 
axial, radial and circumferential velocity 
pressure 

turbulence energy 
energy dissipation 

effective, laminar and turbulent viscosity 

density 

time 

axial, radial and circumferential coordinate 
rectangular coordinate directions 
rotor displacements from its steady-state position 
radial coordinate after transformation 
production term in k-e-model 
constants of the k-e-model 

constants of the k-e-model 

general variable standing for u,v,w,p,k e 
general source term 

nominal seal clearance 

steady-state clearance for an eccentric shaft 

varying seal clearance for orbiting shaft 
radius of the precession motion of the shaft 


perturbation parameter 

rotational frequency of the shaft = RFM*ir/30 

revolutions per minute 

precession frequency of the shaft 

entrance lost-coefficient 

Length of the seal 

stator radius (Fig. 1) 

nominal stator radius (Fig. 1) 

shaft radius (Fig. 1) 

change in the clearance 
eccentricity 


“En 

^Res 


average circumferential velocity 
average axial velocity at the entrance 

sump pressure 
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APRNDIX A: Constants of the k-e. model. 


= 0.09 C t = 1.44 C 2 = 1.92 a k = 1.0 = 1.3 

and G is a given by 
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APPENDIX B: Terms of first order source term. 
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